HERMITE AND BERNSTEIN STYLE BASIS FUNCTIONS FOR CUBIC 
SERENDIPITY SPACES ON SQUARES AND CUBES 



ANDREW GILLETTE 

Abstract. We introduce new Hermite-style and Bernstein- style geometric decompo- 
sitions of the cubic order serendipity finite element spaces S3{P) and 53(/^), as defined 
in the recent work of Arnold and Awanou [Found. Comput. Math. 11 (201 1), 337-344]. 
The serendipity spaces are substantially smaller in dimension than the more conmionly 
used bicubic and tricubic Hermite tensor product spaces - 12 instead of 20 for the square 
and 32 instead of 64 for the cube - yet are still guaranteed to obtain cubic order a pri- 
ori error estimates when used in finite element methods. The basis functions we define 
have a canonical relationship both to the finite element degrees of freedom as well as 
to the geometry of their graphs; this makes the bases ideal for applications employing 
isogeometric analysis where domain geometry and functions supported on the domain 
are described by the same basis functions. Moreover, the basis functions are linear com- 
binations of the commonly used bicubic and tricubic polynomial Bernstein or Hermite 
basis functions, allowing their rapid incorporation into existing finite element codes. 



1. Introduction 

Serendipity spaces offer a rigorous means to reduce the degrees of freedom associated 
to a finite element method while still ensuring optimal order convergence. The 'serendip- 
ity' moniker came from the observation of this phenomenon among finite element prac- 
titioners before its mathematical justification was fully understood; see e.g. [14, 10, 11, 
7]. Recent work by Arnold and Awanou [1,2] classifies serendipity spaces on cubical 
meshes in n > 2 dimensions by giving a simple and precise definition of the space of 
polynomials that must be spanned, as well as a unisolvent set of degrees of freedom for 
them. For degree r convergence, the serendipity space is defined as the span of m mono- 
mials in n variables and is denoted Sr{I^)} Since these m monomials bear no obvious 
correspondence to the m domain points of the n-dimensional unit cube nor to the m 
degrees of freedom described in the papers, alternative local bases for the serendipity 
spaces must be derived. 

In this paper, we provide two coordinate-independent geometric decompositions for 
both 53(7^) and 53(7^), the cubic serendipity spaces in two and three dimensions, re- 
spectively. We present sets of polynomial basis functions associated in a natural and 
symmetric fashion to the cubic serendipity domain points of the square or cube and 
prove that they provide a basis for the corresponding cubic serendipity space. Each basis 
is designated as either Bernstein or Hermite style, as each function restricts to one of 
these common basis function types on each edge of the square or cube. 

The Hermite style bases have the desirable property of being Lagrange-like with re- 
spect to vertex values and coordinate directional derivative values at vertices. More pre- 
cisely, if a Hermite style function is associated to a vertex, it has value 1 at that vertex. 

Date: August 30, 2012. 

Key words and phrases. Serendipity spaces, Hermite and Bernstein functions, isogeometric analysis. 
Supported in part by NSF Award 0715146 and the National Biomedical Computation Resource. 
^The definition of Sr{I^) and its dimension m is reviewed in Section 2.1. 

1 



2 



A. GILLETTE 



value at other vertices, and coordinate directional derivatives equal to at all vertices. 
If a Hermite style function is associated to an edge domain point, it has directional deriv- 
ative equal to 1 at the nearest vertex (in the direction of that domain point) with all other 
coordinate directional derivative values equal to (at vertices) and all vertex function 
values equal to zero. The Bernstein style basis functions have the more simply-stated 
property of obtaining a unique maximal value at the associated domain point. Graphs of 
both types of functions over P are shown in Figures 4 and 5. 




Figure 1 . Cubic serendipity functions on P from [15]. The left function 
is associated to the vertex below the peak. The middle and right functions 
are associated to the edge y = —1 but do not correspond to the domain 
points (±|, — 1) in any canonical or symmetric fashion, making them less 
useful for geometric modeling or isogeometric analysis. 



To the author's knowledge, the only basis functions previously available for cubic or- 
der serendipity finite element purposes employ Legendre polynomials, which lack a clear 
relationship to the domain points. Definitions of these basis functions can be found in 
Szabo and Babuska [15, Sections 6.1 and 13.3]; the two functions from [15] associated to 
the edge ^ = — 1 of are shown in Figure 1 (middle and right). The restriction of these 
functions to the edge gives an even polynomial in one case and an odd polynomial in the 
other, forcing an ad hoc choice of how to associate the functions to the corresponding do- 
main points (±|, — 1). On the other hand, as was just described, the functions presented 
in this paper do have a natural correspondence to the domain points of the geometry. 

Maintaining a concrete and canonical relationship between domain points and basis 
functions is an essential component of the growing field of isogeometric analysis (IGA). 
One of the main goals of IGA is to employ basis functions that can be used both for 
geometry modeling and finite element analysis, exactly as we provide here for cubic 
serendipity spaces. Each function is a linear combination of bicubic or tricubic Bern- 
stein or Hermite polynomials; the specific coefficients of the combination are given in 
the proofs of the theorems. This makes the incorporation of the functions into a vari- 
ety of existing application contexts relatively easy. Note that tensor product bases in 
two and three dimensions are commonly available in finite element software packages 
(e.g. deal.II [4]) and cubic order tensor products in particular are commonly used both in 
modern theory (e.g. isogeometric analysis [9]) and applications (e.g. cardiac electrophys- 
iology models [16]). Hence, a variety of areas of computational science could directly 
employ the new cubic order serendipity basis functions presented here. 

The benefit of serendipity finite element methods is a significant reduction in the com- 
putational effort required for optimal order (in this case, cubic) convergence. Cubic 
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serendipity methods on meshes of squares requires 12 functions per element, an im- 
provement over the 16 functions per element required for bicubic tensor product meth- 
ods. On meshes of cubes, the cubic serendipity method requires 32 functions per element 
instead of the 64 functions per element required for tricubic tensor product methods. Us- 
ing fewer basis functions per element reduces the size of the overall hnear system that 
must be solved, thereby saving computational time and effort. An additional computa- 
tional advantage occurs when the functions presented here are used in an isogeometric 
fashion. The process of converting between computational geometry bases and finite 
element bases is a well-known computational bottleneck in engineering applications [8] 
but is easily avoided when basis functions suited to both purposes are employed. 

The outline of the paper is as follows. In Section 2, we fix notation and summarize 
relevant background on Hermite and barycentric basis functions as well as serendipity 
spaces. In Section 3, we present polynomial Bernstein and Hermite style basis functions 
for S3{P) that agree with the standard bicubics on edges of P and provide a novel geo- 
metric decomposition of the space. In Section 4, we present polynomial Bernstein and 
Hermite style basis functions for Ss{I^) that agree with the standard tricubics on edges 
of reduce to our bases for P on faces of and provide a novel geometric decom- 
position of the space. Finally, we state our conclusions and discuss future directions in 
Section 5. 



2.1. Serendipity Elements. We first review the definition of serendipity spaces and 
their accompanying notation from the work of Arnold and Awanou [1,2]. 

Definition 2.1. The superUnear degree of a monomial in n variables, denoted sldeg(-). 



In words, sldeg(g) is the ordinary degree of q, ignoring variables that enter linearly. 
For instance, the superlinear degree of xy'^z^ is 5. 

Definition 2.2. Define the following spaces of polynomials, each of which is restricted 
to the domain = [-1, 1]^ C W: 

Vr{I^) '= spauj^ {monomials in n variables with degree < r} 

Sr{I^) := spauj^ {monomials in n variables with superUnear degree < r} 

Qr{I^) '= spauj^ {monomials in n variables with each variable degree < r} 

Note that Vr{I^) C Sr{I^) C Qr{I^), with proper containments when r, n > 1. The 
space Sr{I^) is called the degree r serendipity space on the n-dimensional unit cube 
In the notation of the more recent paper by Arnold and Awanou [2], the serendipity spaces 
discussed in this work would be denoted 5^A°(/^), indicating that they are differential 
0-form spaces. The spaces have dimension given by the following formulas (cf. [1]). 



2. Background and Notation 



is given by 





dim Qr{r) 



(r + 1) 



n 
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We write out standard bases for these spaces more precisely in the cubic order cases of 
concern here. 

Vsil^) = span{ 1 , , x\^y\xy , } (2.2) 

linear quadratic cubic 

<S3(/')=P3(/')Uspan{ ^^y^ } (2.3) 

superlinear cubic 

23(/') = S^{I^) U span{x2^2^ xY. . x^V^} (2.4) 
Observe that the dimensions of the three spaces are 10, 12, and 16, respectively. 
7^3(7^) = spanjl, X, y, z, xy, xz, yz, x^ , y^ , z^ , x^y, x^z, xy^ , y^z, xz^ , yz^ , xyz} 



linear quadratic cubic 

(2.5) 



03(7 ) = 7^3(7 ) U spanjx y^x z^y z^xy , ^yz , x yz^ xy x^z , x ^z, x^ x^z } 



superlinear cubic 

(2.6) 

Q3(/') = 53(7^) U span{x^2/^ . . . , x^z'} (2.7) 

Observe that the dimensions of the three spaces are 20, 32, and 64, respectively. 

The serendipity spaces are associated to specific degrees of freedom in the classical 
finite element sense. For a face / of 7^ of dimension > 0, the degrees of freedom 
associated to / for 5^(7^) are (cf. [1]) 



u\ — y uq, qe Vr-2d{f)' 



For the cases considered in this work, n = 2 or 3 and r = 3 so the only non-zero degrees 
of freedom are when / is a vertex {d = 0) or an edge {d = 1). Thus, the degrees of 
freedom for our cases are the values 



u{v), 



J u dt, and j ut dt, (2.8) 



for each vertex v and each edge e of the square or cube. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2. The Bernstein basis and the Hermite basis [i/j] on [0, 1]. 



2.2. Cubic Bernstein and Hermite Bases. For cubic order approximation on square or 
cubical grids, tensor product bases are typically built from one of two alternative bases 
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for7^3([0,l]): 



A. 



1 — 3x + 



■ + 



X 



■ x*" 
x^ 



1 - 3x2 + 
X — 2x2 + x^ 

Jl ^3 



X 

3x2 



- X" 

2x^ 



(l-x)^" 

^1 — x)2x 

^1 — x)x2 
.3 



X 



The set is the cubic Bernstein basis and the set [^] is the cubic Hermite basis. Both 
bases are shown in Figure 2. Bernstein functions have been used recently to provide a 
geometric decomposition of finite element spaces over simplices [3]. Hermite functions, 
while more common in geometric modeling contexts [12] have also been studied in finite 
element contexts for some time [6]. 

The Hermite functions have the following important property relating them to the ge- 
ometry of the graph of their associated interpolant: 



VmGP3([0,1]). (2.9) 



We have chosen these sign and basis ordering conventions so that both bases have the 
same symmetry property: 



/3k{l -x) = k-k{x), ipk{^ -x)= i)^-k{x). 



(2.10) 



The bases and [ip] are related by [13] = V[^] and [ip] = where 



V = 



1-300 

10 

10 

0-31 



13 
10 
10 
3 1 



(2.11) 



Bases for Q3([0, 1]") are easily constructed by taking tensor products of bases for 
Vr{[0, 1]). For instance, note that ® [P] (denoted [/^^j) is a basis for Q3([0, 1]^) where 

[13] ® [/5] = [/3^] = mx)/3,{y) : i,j e {1,2,3,4} } 



The basis function (3i{x)l3j{y) will be denoted I3ij. Similarly, [13] ® [0] ® [0] (denoted 
[/3^]) is a basis for Q3([0, 1]^) with basis functions denoted /5jjfc. 

We now rephrase the fact that [0] is a basis for 7^3 ([0, 1]) as a 'precision property,' i.e. 
that linear combinations of the basis functions recover certain monomials exactly. The 
extension of this property to tensor product bases will be used to prove that the sets of 
functions presented in Sections 3 and 4 are indeed bases for S2,{P) and <S3(/^). 
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Proposition 2.3. For < s^t < 3, the precision properties of [/3], and [/3^] take 
on the respective forms 

4 



r 

X 



E (!:>.. (2.12) 



4 4 



EEE(r-;)C:;)(::>* 



Prc^o/ One easily confirms that = ^04, x'^ = f^^ + (3^, x^ = f^^ + 2(3^ + /32, and 
= /54 + 3;53 + 3/32 + /3i, which is the statement of (2.12). The other two statements 
follow immediately from the first. □ 

We have a similar proposition for the Hermite basis [ip] and its tensor products [ijj'^] := 
[ip] (8) [ip] and [i/j^] := [^0] [ip] ['0]. As before, ipij means ipi{x)ipj{y) and ipijk means 
^i{x)ijj{y)^k{z). 

Proposition 2.4. Let 

^^.-=E(4I1)^- (2.15) 

where Vai denotes the (a, i) ^fzfry frc^w, column) of N from (2.11). For < r, 5, t < 3, the 
precision properties oflip], [ip'^], and [ip^] take on the respective forms 

4 



r 

X 



^Sr^i^pi, (2.16) 



i=l 
4 4 



^"^Sr^iSsJ^pij, (2.17) 

4 4 4 

i=i j=i k=i 

Proof. By swapping the order of summation, we see that 



E ^r,.i^^ = E (4 _ I) (e = E (L 1 



by (2.12), proving (2.16). The other two statements follow immediately from the first. 

□ 

Transforming the bases [f3] and [ip] to domains other than [0, 1] is straightforward. If 
T : [a, b] [0, 1] is linear, then replacing x with T{x) in each basis function expression 
for [/3] and [ip] gives bases for V3{[a,b]). Note, however, that the derivative interpolation 
property for [ip] must be adjusted to account for the scaling: 

u = u{l)iJi{T{x)) + \b- a\u\0)ij2{T{x)) 

-\b- a\u\l)iJs{T{x)) + u{l)^4{T{x)), G V^{[a, b]). 

(2.19) 
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In geometric modeling applications, the coefficient |6 — a| is sometimes left as an ad- 
justable parameter, usually denoted s for scale factor. [5, 16]. For all the Hermite and 
Hermite style functions, we will use derivative-preserving scaling which will include 
scale factors on those functions related to derivatives; this will be made explicit in the 
various contexts where it is relevant. 

Remark 2.5. Both [/3] and [ip] are Lagrange like at the endpoints of [0, 1], i.e. at an 
endpoint, the only basis function with non-zero value is the function associated to that 
endpoint (/3i or ipi for 0, /34 or 7/^4 for 1). This means the two remaining basis functions of 
each type (/52, (3s or ?/;2, ^3) are naturally associated to the two edge degrees of freedom 
(2.8). We will refer to these associations between basis functions and geometrical objects 
as the standard geometrical decompositions of and [i/j]. 



3. Local Bases for Ss{P 





Figure 3. On the left, ordered pair indices from X are shown next to 
the domain point of [0, 1]^ to which they correspond. On the right, only 
those ordered pairs used for the serendipity basis are shown. The corre- 
spondences V ^ vertices, E ^ edge points, and D ^ domain interior 
points are evident. 



Before defining local bases on the square, we fix notation for the domain points to 
which they are associated. For [0, 1]^, define the ordered paired index set 

X:= {{t,j}\i,je{l,...,4}}- 
Then X is the disjoint union V U E U D where 



V 


■■={{hJ 


}\tje{lA}}; 


(3.1) 


E 


■■={{hJ 


} exactly one of {i, j} is in {1, 4}} ; 


(3.2) 


D 


■■={{hJ 


}\t,je{2,3}} 


(3.3) 



The V indices are associated with vertices of [0, 1]^, the E indices to edges of [0, 1]^, and 
the D vertices to the domain interior to [0, 1]^. The relation between indices and domain 
points of the square is shown in Figure 3. 



3.1. A local Bernstein style basis for Ss{P). Under the notation and conventions es- 
tablished in Section 2, we establish a local Bernstein style basis for S3{P) where / := 
[—1,1]. Define the following set of 12 functions, indexed by V U E; note that each 
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function is to be scaled by a factor of 1/16. 

{x -l)iy- l)(-2 -2x + x'^ -2y + y^) 
-{x - l){y + l)(-2 -2x + x^ + 2y + y"^) 
-{x + l){y - l)(-2 + 2x + - 2|/ + y"^) 
{x + l)(y + l)(-2 + 2a; + + 2y + y'^) 
-(x-l)(^-l)2(y + l) 
{x-l){y-l){y + lf 
(x + l)(^-l)2(y + l) 
-{x + l){y-l){y + lf 
-{x-lf{x + l){y-l) 
{x-l){x + lf{y-l) 
{x-lf{x + l){y + l) 
-{x-l){x^l)\y + l) 

Fix the basis orderings 

[?^] •= [6ii 64,61164, 62,63,62,63,61,61,64,64], 



























S42 




63 




61 




61 








_<^34_ 





11 



indices in V 

) Pu^ ^41, 



^44, /3: 



12 



indices in E 
^13, /342, ^43^ ^21, ^31, /324, /534, 



1 

16 



(3.4) 



(3.5) 



/322,/323,^32,^33] (3.6) 



indices in V 



indices in E 



indices in D 




^11 ^21 <^31 



Figure 4. The top row shows 3 of the 16 bicubic Bernstein functions on 
P while the bottom row shows 3 of the 12 cubic Bernstein style serendip- 
ity functions. The visual differences are subtle, although some changes in 
concavity can be observed. Note that functions in the same column have 
the same values on the edges of P. 

Theorem 3.1. Let /3}^ denote the scaling of to i.e. 

/3L:=A((a; + l)/2)/3^((?/ + l)/2). 
The set [^^] has the following properties: 
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(i) [.^2] isabasisfor<S3(/2). 

(ii) For any £m eVuE, ^^rn is identical to on the edges of P. 

(iii) [^^] is a geometric decomposition of S'i{P). 

Proof. For (i), we scale [^'^] to [0, 1]^ to take advantage of a simple characterization of the 
precision properties. Let [^2j[o,i] denote the set of scaled basis functions ^f^\x, y) := 
iimi'^x — 1, 2?/ — 1). Given the basis orderings in (3.5)-(3.6), it can be confirmed directly 
that [C2][°'il is related to [jS"^] by 

where B is the 12x16 matrix with the structure 



B := [ I I B' ] , 

where I is the 12 x 12 identity matrix and B' is the 12 x 4 matrix 



(3.8) 



-4 -2 -2 -1" 

-2 -4 -1 -2 

-2 -1 -4 -2 

-1 -2 -2 -4 



2 

1 

2 

1 





2 

1 
1 

2 




1 

2 


2 

1 




1 

2 

1 

2 



(3.9) 



Using iji G X to denote an index for and £m G V U £' to denote an index for 
the entries of B can be denoted by 6f J* so that 



[0,1] 
tm ' 





^11 


. . • 


• /)11 1 
^33 




llm 
^11 


iJtm 
ij 


' ^33 




_ ^11 




^33 _ 



(3.10) 



We now observe that for each ij e X, 



3 — r 
4-i 



3- s 

4- J 



3 — r 
4-£ 



3-s 



(3.11) 



for all (r, 5) pairs such that sldeg(x^^^) < 3 (recall Definition 2.1). Note that this claim 
holds trivially for the first 12 columns of B, i.e. for those ij ^ V E d X. For 
ij G D C X, (3.1 1) defines an invertible linear system of 12 equations with 12 unknowns 
whose solution is the ij column of B'; the 12 (r, s) pairs correspond to the exponents of x 
and y in the basis ordering of Ss{P) given in (2.2)-(2.3). Substituting (3.11) into (2.13) 
yields: 



E E 

ijex \imevuE 




10 



A. GILLETTE 



Swapping the order of summation and regrouping yields 



emeVuE ^ ^ ^ ^ \ijex 
The inner summation is exactly ^f^^ by (3.7), implying that 

^3 - r\ / 3 

A-e U 



= E 



imeVUE 



m 



cm 



(3.12) 



for all (r, s) pairs with sldeg(x^^^) < 3. Since has 12 elements which span the 

12 dimensional space 53([0, 1]^), it is a basis for 53([0, 1]^). By scaling, [^^] is a basis 

For (ii), note that an edge of [0, 1]^ is described by an equation of the form {x or y} = 
{0 or 1}. Since /52(t) and /33(t) are equal to at t = and t = 1, = on the edges of 
[0, 1]^ for any ij G D. By the structure of B from (3.8), we see that for any ^m G ^ U 



^[0.1] 



Thus, on the edges of [0, 1]^, and are identical. After scaling back, we have 
and identical on the edges of as desired. 

For (iii), the geometric decomposition is given by the indices of the basis functions, 
i.e. the function (^^^ is associated to the domain point for Im ^ V [J E. This follows 
immediately from (ii), the fact that is a tensor product basis, and Remark 2.5. 

□ 

Remark 3.2. It is worth noting that the basis [^^] was derived by essentially the reverse 
order of the proof of part (i) of the theorem. More precisely, the twelve coefficients in 
each column of B define an invertible linear system given by (3.1 1). After solving for the 
coefficients, we can immediately derive the basis functions via (3.7). By the nature of this 
approach, the edge agreement property (ii) is guaranteed by the symmetry properties of 
the basis [(5] . This technique was inspired by a previous work for Lagrange-like quadratic 
serendipity elements on convex polygons [13]. 

3.2. A local Hermite style basis for Sz{P). We now establish a local Hermite style 
basis [^^] for S'i{P) using the bicubic Hermite basis [?/;^] for Q3([0, 1]^). Define the 
following set of 12 functions, indexed by 1/ U note that each function is to be scaled 
by a factor of 1/8. 

-(x - l){y - l)(-2 + x + x'^ + y + y^y 
[x - l){y + l)(-2 + X + - 2/ + 2/2) 
{x + l){y - l)(-2 - X + + 2/ + 2/2) 
-(x + l)\y + l)(-2 -x + x^ -y + y^) 
-(x -!)(?/- 1)2(2/ + 1) 
(x- !)(?/- !)(?/ + 1)2 1 
(x + !)(?/- 1)2(^ + 1) 
-{x + l){y-l){y + lf 

— (x — l)2(x 

(x — l)(x 4 

(x — l)2(x 

— (x — l)(x 



ijeD 



(3.13) 











^41 




^44 








As 




^2 




^3 




^21 




^31 




^24 




^34 





i){y-i) 
my -I) 



(3.14) 
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Fix the basis orderings 

[^'] := [^11,^14,^41,^44, ^12,^13,^42,^43,^21,^31,^24,^34], (3.15) 

^ V ' ^ V ' 

indices in V indices in E 

[^2] [^^^,^^4,^4^,^44, ^12, ^13, ^42, ^43, ^21, ^31, ^24,^^34, ^22,^23,^32,^33] 

" V ' " V ' " V ' 

indices in V indices in E indices in D 

(3.16) 




Figure 5. The top row shows 3 of the 16 bicubic Hermite functions on 
1^ while the bottom row shows 3 of the 12 cubic Hermite style serendipity 
functions. The visual differences are subtle, although some changes in 
concavity can be observed. Note that functions in the same column have 
the same values on the edges of 



Theorem 3.3. Let "^{^ denote the derivative-preserving scaling of "^^^ to i.e. 

- + l)/2)^^((2/ + l)/2), im G V, 

:= 2^,{{x + \)l2)^^{{y + l)/2), im G E, 

The set \d^\ has the following properties: 

(i) \^^\ isabasisfor53(/2). 

(ii) For any im G 1/ U is identical to "^{^ on the edges of 

(iii) \d^\ is a geometric decomposition of S'^{T^\ 

Proof. The proof follows that of Theorem 3.1 so we abbreviate proof details that are 
similar. For (i), let [i?^]!^'^] denote the derivative-preserving scaling of \&^\ to [0, 1]^; the 
scale factor is 1/2 for functions with indices in E. Given the basis orderings in (3.15)- 
(3.16), we have 

=]H[[7/;2] (3.17) 
where H is the 12 x 16 matrix with the structure 

M := [ I I ] , (3.18) 



12 
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where I is the 12 x 12 identity matrix and is the 12 x 4 matrix 



■ -1 


1 


1 


-1 ■ 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


-1 


1 


1 


-1 


-1 





1 








-1 





1 


1 





-1 








1 





-1 


-1 


1 














-1 


1 


1 


-1 














1 


-1 



(3.19) 



Denote the entries of H by hfj^ (cf. (3.10)). Recalling (2.15), observe that for each 

ij e X, 

^r^i^sj ^ ^ ^r,i^s,jh— , (3.20) 

for all (r, 5) pairs such that sldeg(x^j/^) < 3. Similar to the Bernstein case, we substitute 
(3.20) into (2.17), swap the order of summation and regroup, yielding 



imeVUE 



The inner summation is exactly i^f^^ by (3.17), implying that 



[0,1] 



(3.21) 



imeVUE 



for all (r, s) pairs with sldeg(x^^^) < 3, proving that [i?^]!^'^] is a basis for Ss{[0, 1]^). By 
derivative-preserving scahng, is a basis for S3{P). 

For (ii), observe that for any ij G ipij = on the edges of [0, 1]^ by virtue of the 
bicubic Hermite basis functions' definition. By the structure of H from (3.18), we see 
that for any Im^V \JE, 



[0,1] 



ijeD 



(3.22) 



Thus, on the edges of [0, 1]^, "df^^ and ?/;^^ are identical. After scaling back, we have 
i^irn and ijjj^ identical on the edges of as desired. 

For (iii), the geometric decomposition is given by the indices of the basis functions, 
i.e. the function i^i^ is associated to the domain point for £m e V U E. This follows 
immediately from (ii), the fact that [ip'^] is a tensor product basis, and Remark 2.5. □ 

4. Local Bases for Ss{I^) 

Before defining local bases on the cube, we fix notation for the domain points to which 
they are associated. For [0, 1]^, define the ordered triple index set 



Y :={{i,j,k}\i,j,ke{l,...A}}- 



HERMITE AND BERNSTEIN STYLE BASIS FUNCTIONS FOR CUBIC SERENDIPITY SR\CES 13 



244 344 444 244 344 444 



144 
143 


\^ 234 
134\^* 

124\ 


334\444 
224 324\ 


444 




144 
143 


134 


-• — 

124 


— •- 






s444 


444 


142 


133 114 
• 

123 

132* 113 


214 314 

• • 

213 313 


414 
413 




142 






114 
113 




214 


314 


414 
413 ^ 


141 


\^122* 

13l\^112 


• • 

212 312 


412 




141 


131 




112 








412 ^ 




12l\ 












121 














111 


i • •- 

211 311 


411 










111 




•- 

211 


•- 

311 


411 



Figure 6. On the left, ordered triple indices from Y are shown next to 
the domain point of [0, 1]^ to which they correspond. Points hidden by 
the perspective are not shown. The origin is at the point labelled 111; 
the positive x, y and z axes go right, back, and up, respectively. On the 
right, only those indices used for the serendipity basis are shown. The 
correspondences V ^ vertices, E ^ edge points, F ^ face interior 
points, and M ^ domain interior points are evident. 



Then Y is the disjoint union V U E U F U M where 



V:={{i,j,k}\i,j,ke{lA}}'. (4.1) 

E := {{ij, k} I exactly two of {ij, k} are in {1, 4}} ; (4.2) 

F := {{i, j, k} I exactly one of {i, j, k} is in {1, 4}} ; (4.3) 

M:={{i,j,k}\i,j,ke{2,3}} (4.4) 



The V indices are associated with vertices of [0, 1]^, the E indices to edges of [0, 1]^, the 
F indices to face interior points of [0, 1]^, and the M vertices to the domain interior of 
[0,1]^. The relation between indices and domain points of the cube is shown in Figure 6. 



4.1. A local Bernstein style basis for Ss{I^). Under the notation and conventions es- 
tablished in Section 2, we are ready to establish a local Bernstein style basis for S3{P) 
where / := [—1,1]. Define the following set of 32 functions, indexed hy V U E C Y; 



14 



A. GILLETTE 



note that each function is to be scaled by a factor of 1/32. 
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S413 
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-{x + 1 


<C424 
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<C431 
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-{x + l 


<^442 
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-{x + l 



-2x + x^ -2y + y^ -2z + z^)' 
-5 - 2x + x^ - 2y + y^ + 2z + z'^) 
-5 - 2x + x^ + 2y + y^ - 2z + z^) 
(-5 -2x + x^ + 2y + y^ + 2z + z^) 
-5 + 2x + x'^ -2y + y'^ -2z + z^) 
(-5 + 2x + x'^ -2y + y^ + 2z + z'^) 
(-5 + 2a; + + 2|/ + ?/2 _ 2z + z'^) 
-5 + 2x + x'^ + 2y + y'^ + 2z + z'^) 
y-l)iz-inz + l) 
{y - l){z - l){z + If 
y-iny+l){z-l) 

(y-lf{y+l){z + l) 

iy-l){y + inz-l) 

y-l){y+lf{z+l) 

{y + l){z-inz + l) 
y + l)iz-l)iz + lf 

ix + l)iy-l)iz-l) 
\x + l){y-l){z + l) 
^{x + l){y + l){z-l) 
{x + l){y + l){z + l) 
{x + iny-l)iz-l) 
X + l)\y - l){z + 1) 

x + lf{y + l){z-l) 
{x + lf{y + l){z + l) 
(y-l)(z-inz + l) 

y-l)(z-l){z+ir 

{y-iny + l){z-l) 
y-iny+l){z + l) 
y-l){y + lf{z-l) 
{y-l){y + lfiz + l) 
y + l){z-lf{z + l) 
{y + l){z-l){z + ir 



1 

32 



(4.5) 



We fix the following basis orderings, with omitted basis functions ordered lexicographi- 
cally by index. 



111) 



444, <Cll2, 



443 



indices in V indices in E 

[ All, . . . , /5444, /5ll2, • • • , /5443, /5l22, • • • , ^433^ ^222^ • • • , /^333^ 



(4.6) 
(4.7) 



indices in V 



indices in E 



indices in F 



indices in M 



Theorem 4.1. Let /?/^^ denote the scaUng of Pimn to i.e. 

^Ln := + l)/2)^m((2/ + l)/2)M{z + l)/2). 
The set [^^] has the following properties: 

(i) [^^] is a basis for Ss{P). 

(ii) [^^] reduces to [(^^] on faces of P. 

(iii) For any ^mn eVUE, ^£rnn is identical to on edges of 
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(iv) [^^] is a geometric decomposition of Ss{I^). 

Proof. The proof is similar to that of Theorem 3.1 so we abbreviate some details. For (i), 
let [^^][°'il denote the scaling of [<^^] to [0, 1]^. Given the basis orderings in (4.6)-(4.7), 
we have 

[^3] [0,1] =u[/3'] (4.8) 

where U is the 32 x 64 matrix with the structure 

U := [ I I ] , (4.9) 

where I is the 32 x 32 identity matrix and is a specific 32 x 32 matrix whose entries 
are integers ranging from -16 to 4. Instead of writing out in its entirety, we describe 
its properties and how it can be constructed (cf. Remark 3.2). 

Using ijk to denote an index for f3ijk and £mn eVUEcYio denote an index 
for ^f^l, the entries of U will be denoted by kf^^ so that 





^111 


• 7/111 . . 


• 7/111 1 
^333 




^111 




^333 




_ ^111 




^333 J 



(4.10) 



The columns of U satisfy the relationship 




(4.11) 



for all (r, s, t) tuples such that sldeg(x^7/^z^) < 3. This property defines the entries of U 
uniquely since for each ijk eY it gives an invertible linear system of 32 equations with 
32 unknowns whose solution is the ijk column of U. The (r, s, t) tuples should be taken 
in the order given in (2.5)-(2.6). See Remark 3.2 and the text after (3.11) for more on 
this process. 

As in previous proofs, regrouping and recognizing an expression for ^f^l gives 




(4.12) 



proving, after scaling, that [^^] is a basis for S3{P). 

For (ii), the claim can be confirmed directly by calculation, e.g. ^iii(x, t/, — 1) = 
y) or ^i42(^, 1, z) = Ci2(x, z), etc. 

For (iii), note that an edge of [0, l]^ is described by two equations of the form {x, y, or z} = 
{0 or 1} where two distinct variables must be chosen for the two equations. Since /32{t) 
and /3s{t) are equal to at t = and t = 1, /3ijk = on the edges of [0, 1]^ for any 
ijk G M. Further, for ijk G F, without loss of generality, assume that i G {1,4} so 
that j, k G {2, 3}. Since every edge is described by at least one equation of the form 
{y or z} = {0 or 1}, either f3j{y) or f3k{z) is identically zero on every edge. Thus, for 
ijk G F U M, f3ijk = on the edges of [0, l]^. 

By the structure of U from (4.9), we see that for any £mn gVUE, 

a = ^imn+ E '^^^^>^- (4.13) 

ijkeFUM 
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Thus, on the edges of [0, 1]^, ^f^l and Pimn are identical. After scahng back, we have 
^irnn ^rid identical on the edges of as desired. 

For (iv), the geometric decomposition is given by the indices of the basis functions, 
i.e. the function (^^^^ is associated to the domain point for £mn ^ V U E. This follows 
immediately from (ii) and (iii), the fact that is a tensor product basis, and Remark 2.5. 

□ 



4.2. A local Hermite style basis for Ss{P). We now estabhsh a local Hermite style 
basis [d^] for Ss{I^) using the tricubic Hermite basis [ip^] for Q3([0, 1]^). Define the 
following set of 32 functions, indexed hy V U E; note that each function is to be scaled 
by a factor of 1/16. 



'^111 








^Ul 




^144 




^411 








1^441 








^112 




^113 




^121 




^?124 




^131 




1^134 




1^142 




^143 




^211 




'^214 




'^241 




'&244 




^311 




1^314 




1^341 




1^344 




1^412 








^?421 




^?424 




^?431 




^^434 




1^442 




.^443. 





{x - l){y - l)iz - 1) 
-ix-l)iy-l)iz + l 
-ix-l)iy + l)iz-l 
{x-l){y + l){z + l) 
-{x + l){y-l){z-l 
{x + l){y-l){z + l) 
{x + l){y + l){z-l) 
-{x + l){y + l){z + l 
{x-1 

— {x — . 
{x-1 

— {x — . 

— {x — '. 
{x - 

— {x 
{x - 
{x-l 

— ix — 



-2 + X + + ?/ + ?/2 + z + 2;^) ' 
{-2 + X + x"^ ^-y + y'^ - z + z"^) 
{-2 + x + x'^ -y^-y"^ + z + z"^) 
—2 + x + x'^ — y + y'^ — z + z'^) 
- x + x'^ + y + y'^ + z + z'^) 
X + x'^ + y + y'^ — z + z'^) 
X + x'^ — y + y'^ + z + z'^) 



(-2 

-2- 

-2- 



1 



1 



— (a; — 
{x-1 

— {x — 
{x-1 
{x-1 

— {x — 
-{x + 

{x + 1 
-{x + 
{x + 1 

-{x + 1 
{x + 1) 
-{x + 1 



-2 — X + x'^ — y + 
l){z-lf{z + l) 
-l){z-l){z + lf 
lf{y + l){z-l) 

-lf{y + l){z + l) 

-l){y + lf{z-l) 
l){y + lf{z + l) 
Vl){z-lf{z + l) 



z + z^) 



y - 

{y 
y- 
{y 
(y 
y - 
iy 

y + l){z-l){z + lf 
{x + l){y-l){z-l) 
'{x + l){y-l){z + l) 
'{x + l){y + l){z-l) 
{x + l){y + l){z + l) 
{x + lf{y-l){z-l) 
x + lf{y-l){z + l) 
x + lf{y + l){z-l) 
{x + lf{y + l){z + l) 
{y-l){z-inz + l) 
y-l){z-l){z + ir 

{y-lf{y + l){z-l) 

y-lf{y + l){z + l) 
y-l){y + lf{z-l) 

iy-i){y 

y + l){z- 
{y + l){z 



■mz + l) 

mz+i) 

l){z + iy 



1 

16 



(4.14) 
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We fix the following basis orderings, with omitted basis functions ordered lexico- 
graphically by index. 

[l9^]:=[l9ni,...,^444, ^112,...,^443 ], (4.15) 
indices in V indices in E 

:= . . . ,?/;444, ^0112, • • • , ^443, ^0122, . . . ,^433, fe,...,^333, ] (4.16) 

^ V ' ^ V ' ^ V ' ^ V ' 

indices in V indices in E indices in F indices in M 

Theorem 4.2. Let ^/^^ denote the derivative-preserving scaling of i/jimn to i.e. 

:= Mi^ + l)/2)^m((?/ + l)/2)M{^ + l)/2), £mn G V, 

^\^^ := 2^,{[x + \)l2)^m[ky + l)/2)^n((^ + l)/2), imn G E. 

The set has the following properties: 

(i) \&^\ is a basis for 53(/^). 

(ii) reduces to \&^\ on faces of 

(iii) For any imn G 1/ U "dtmn is identical to "^{^^ on edges of 

(iv) \&^\ is a geometric decomposition of S'^{1^\ 

Proof, The proof is similar to that of Theorem 3.3 so we abbreviate some details. For (i), 
let denote the scaling of [d^] to [0, 1]^. Given the basis orderings in (4.15)-(4.16), 

we have 

[^3|[o,i] = (4.17) 

where W is the 32 x 64 matrix with the structure 

W := [ I I W ] , (4.18) 

where I is the 32 x 32 identity matrix and W is a specific 32 x 32 matrix whose entries 
are in {—1, 0, 1}. The matrix W is constructed similarly to the matrix U from the proof 
of Theorem 4.1; the columns of W satisfy the relationship 

^r,i^sj^t,k — ^ ^ ^r,i^sj^t,k^ijk^ ) (4.19) 
imneVUE 

for all (r, s, t) tuples such that sldeg(x^y^z^) < 3. Similar to previous proofs, this yields 

X^'y'z^ = ^ Sr/£s,m£t,n^imn, (4.20) 
ImneVUE 

proving, after derivative-preserving scaling, that [d^] is a basis for Ss{I^). 

For (ii)-(iv), the proof is similar to the corresponding parts of the proof of Theorem 4. 1 . 

□ 

5. CONLUCSIONS AND FUTURE DIRECTIONS 

The basis functions presented in this work are well-suited for use in both geometric 
modeling and finite element applications, as discussed in the introduction. Moreover, 
the proof techniques used for the theorems suggest a number of promising extensions. 
Similar techniques should be able to produce Bernstein style bases for higher polynomial 
order serendipity spaces, although the introduction of interior degrees of freedom that 
occurs when r > 3 requires some additional care to resolve. Some higher order Hermite 
style bases may also be available, although the association of directional derivative values 
to vertices is somewhat unique to the r = 3 case. Pre-conditioners for finite element 
methods employing our bases are still needed. The fact that all the functions we defined 
are fixed linear combinations of standard bicubic or tricubic basis functions suggests that 
such pre-conditioners will have a straightforward construction. 



18 



A. GILLETTE 



References 

[1] D. Arnold and G. Awanou. The serendipity family of finite elements. Foundations of Computational 

Mathematics, ll(3):337-344, 2011. 
[2] D. Arnold and G. Awanou. Finite element differential forms on cubical meshes. Arxiv preprint 

arXiv: 1204.2595, 2012. 

[3] D. Arnold, R. Falk, and R. Winther. Geometric decompositions and local bases for spaces of fi- 
nite element differential forms. Computer Methods in Applied Mechanics and Engineering, 198(21- 
26): 1660-1672, 2009. 

[4] W. Bangerth, R. Hartmann, and G. Kanschat. deal, ii — a general-purpose object-oriented finite ele- 
ment library. ACM Transactions on Mathematical Software (TOMS), 33(4):24-es, 2007. 
[5] C. Bradley, A. PuUan, and R Hunter. Geometric modeling of the human torso using cubic hermite 

elements. Annals of Biomedical Engineering, 25:96-111, 1997. 
[6] P. Ciarlet and P. Raviart. General Lagrange and Hermite interpolation in with applications to finite 

element methods. Archive for Rational Mechanics and Analysis, 46(3): 177-199, 1972. 
[7] P. G. Ciarlet. The Finite Element Method for Elliptic Problems, volume 40 of Classics in Applied 

Mathematics. SIAM, Philadelphia, PA, second edition, 2002. 
[8] J. Cottrell, T. Hughes, and Y. Bazilevs. Isogeometric analysis: toward integration of CAD and FEA. 

John Wiley & Sons Inc, 2009. 
[9] J. Evans and T. Hughes. Explicit trace inequalities for isogeometric analysis and parametric hexahe- 
dral finite elements. Numerische Mathematik, pages 1-32, 2011. 
[10] T. Hughes. The finite element method. Prentice Hall Inc., Englewood Cliffs, NJ, 1987. 
[11] J. Mandel. Iterative solvers by substmcturing for the p- version finite element method. Computer 

Methods in Applied Mechanics and Engineering, 80(1-3): 117-128, 1990. 
[12] M. Mortenson. Geometric Modeling. John Wiley and Sons, 3rd edition, 2006. 
[13] A. Rand, A. Gillette, and C. Bajaj. Quadratic serendipity finite elements on polygons using general- 
ized barycentric coordinates. arXiv: 1109.3259, 2011. 
[14] G. Strang and G. J. Fix. An analysis of the finite element method. Prentice-Hall Inc., Englewood 
Cliffs, N. J., 1973. 

[15] B. Szabo and I. Babuska. Finite element analysis. Wiley-Interscience, 1991. 

[16] Y. Zhang, X. Liang, J. Ma, Y. Jing, M. J. Gonzales, C. Villongco, A. Krishnamurthy, L. R. Frank, 
V. Nigam, P. Stark, S. M. Narayan, and A. D. McCuUoch. An atlas-based geometry pipeline for 
cardiac Hermite model construction and diffusion tensor reorientation. Medical Image Analysis, (in 
press), 2012. 
E-mail address: akgillette@ucsd . edu 



Department of Mathematics, University of California San Diego, La Jolla CA 92093 



